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Domain walls and droplet-like excitation of the random-field Ising magnet are studied in d = 
{3, 4, 5, 6, 7} dimensions by means of exact numerical ground-state calculations. They are obtained 
using the established mapping to the graph-theoretical maximum-flow problem. This allows to study 
large system sizes of more than five million spins in exact thermal equilibrium. All simulations are 
carried out at the critical point for the strength h of the random fields , h = h c (d), respectively. Using 
finite-size scaling, energetic and geometric properties like stiffness exponents and fractal dimensions 
are calculated. Using these results, we test (hyper) scaling relations, which seem to be fulfilled below 
the upper critical dimension d u — 6. Also, for d < d u , the stiffness exponent can be obtained from 
the scaling of the ground-state energy. 
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I. INTRODUCTION 

The random-field Ising magnet (RFIM) is one of 
the most frequently studied models for magnetic sys- 
tems with quenched disorder. For d = 3 and 
higher dimensions^ it is known to undergo a phase 
transition 2 ^— at a critical temperature T c which depends 
on the disorder strength. For low temperatures and weak 
disorder the ferromagnetic interactions dominate and the 
system is ferromagnetically long-range ordered. For large 
temperature or strong disorder, the RFIM exhibits no 
long-range order and behaves like a paramagnet in a field. 

Numerically, the nature of this phase transition can be 
studied in many cases conveniently by means of exact 
ground-state (GS) calculations (see below), based on an 
established mapping to the graph-theoretical maximum- 
flow problem. This allows to treat large system sizes in 
thermal equilibrium, in contrast to Monte Carlo simu- 
lations. For Gaussian disorder, the phase transition is 
of second order along the full transition line and can be 
characterized by critical exponents in the usual way, like 
v, describing the divergence of the correlation length, and 
a, describing the behavior of the specific heat. Never- 
theless, the RFIM behaves differently compared to the 
standard ferromagnet. In particular, the hyper-scaling 
relation dv = 2 — a has to be changed^ by including a 
positive parameter 8, yielding 

v{d-6) = 2-a. (1) 

In general, the value of 9 can be obtained directly from 
the scaling of domain- wall energies, e.g., induced by 
changing the boundary conditions, and is known as stiff- 
ness exponent. The study of such domain-wall exci- 
tations was pioneered in the field of spin glasses i 15 ' 16 
Here, a comprehensive understanding of the nature of 
the behavior of two-dimensional spin glasses could be 
obtained) 17 ' 18 which turned out to be compatible with 
the droplet picture . 15 ' 19 " — 

According the droplet picture, the energy scaling of 
droplet-like excitations should be the same as for do- 



main walls. For two-dimensional spin glasses, this was re- 
cently confirmed via using modified GS algorithms £2r21 
Later on, the value of 9 was also determined in higher 
dimensions^ 5 - - — up to the upper critical dimension, in 
this case via (heuristic) ground-state calculations before 
and after changing the boundary conditions, respectively. 

For the RFIM, which is the subject of this work, 
domain-wall studies similar to the spin-glass case, i.e. 
based on GS calculations, were performed for three- and 
four dimensions* 7 -*^ but according our knowledge not in 
higher dimensions. Droplet type low- or lowest-energy 
excitations were only obtained in three dimensions so 
far j 12 ' 30 For three dimensions at finite temperature, also 
free-energy barrriers where calculated recently 31 ' 32 

Furthermore, the prediction for the RFIM upper crit- 
ical dimension d u — 6 was confirmed 13 via exhaustive 
exact GS calculation up to d = 7. This shows that the 
RFIM can be investigated conveniently by numerical ex- 
act algorithms even close to and above the upper critical 
dimension. 

Hence, it is the purpose of this work to study domain- 
wall and droplet excitations of the RFIM in dimensions 
d = 5,6,7 (and for d = 3,4 for comparison), similar 
to the corresponding d = 2 studies for spin glasses, 
with the striking difference that for the RFIM an exact 
polynomial-time GS algorithm is available for any dimen- 
sion, allowing to treat much larger system sizes of more 
than five million spins in thermal equilibrium. Another 
difference is that we performed the study for the RFIM 
right at the zero-temperature disorder-critical point (like 
the previous work in lower dimensions), since the scaling 
in the ferromagnetic and paramagnetic phases should be 
trivial. We analyzed energetic and geometric properties 
of the excitations using finite-size scaling. We compare 
the results of different excitations, which should agree 
according the droplet picture, and verify the above men- 
tioned hyper-scaling relation. Also we compare the geo- 
metric (fractal) properties of these excitations. Treating 
system sizes up to and above the upper critical dimension 
allows us to observe the transition to mean-field behav- 
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To state the model in detail, the RFIM consists of N 
Ising spins Si — ±1 located on the sites of a hyper-cubic 
lattice with periodic boundary conditions (PBC) in all 
directions. The spins couple to each other and to local 
net fields. Its Hamiltonian reads 



H = -J^2SiSj-^2{hei)Si 



(2) 



It has two contributions. The first covers the spin-spin 
interaction, where J is the ferromagnetic coupling con- 
stant between two adjacent spins and (i, j) denotes pairs 
of next-neighbored spins. The second part of the Hamil- 
tonian describes the coupling to local fields hi = hsi, The 
factor h is the disorder strength and £i the quenched dis- 
order, i.e. Gaussian distributed with zero mean and unit 
width. 

The paper is organized as follows: In section [TT1 we de- 
scribe in principle how the GSs are calculated. The next 
section covers the definition and use of the different exci- 
tations and their theoretical background. Then we state 
our results in section IIVI and finish with our conclusions 
and discussion. 



II. GROUND STATES 

The phase space of the RFIM consists of a ferromag- 
netic and a paramagnetic phase. The transition from 
one phase to the other takes place at a critical point 
Pc = (h c ,T c ). The transition can be triggered vary- 
ing the temperature T or varying the standard devia- 
tion of the disorder distribution h. From Ref. [3^, it is 
known via renormalization group calculations, that the 
critical behavior of the RFIM is controlled from the zero- 
temperature fixed point. Hence, it is possible to focus on 
T = = const, and vary h to study the phase transi- 
tion. Here, we concentrated on T = 0, h f» h c , to study 
excitations right at at the critical point. 

At T = it is possible to calculate exact ground states 
in a very efficient way. Following an approach from Ref. 
I34ll35l . a d-dimensional hypercubic realization of the dis- 
order {hei} can be mapped to a graph with N + 2 nodes 
and (d + 2)N + 1 edges with suitable edge capacities, 
where N is the number of spins of the RFIM. On this 
graph a sophisticated maximum flow/minimum cut al- 
gorithm can be applied. 3 -^ The resulting minimum cut 
directly correspond to the GS spin configuration {Si} of 
that specific realization of the disorder. For our simula- 
tions the implementation of the maximum flow algorithm 
from the LEDA library 3 ^ is used. For the RFIM, the 
actual runtime of the algorithm increases only slightly 
stronger than linear with the number N of spins. 39 



III. DOMAIN WALLS AND DROPLET 
EXCITATIONS 

We studied two types of excitations, domain walls and 
droplets. 

The domain walls treated in this work separate spin 
regions which are effected from changed (boundary) con- 
ditions from unaffected spins. Following Ref. [7| we 
forced boundary spins along distinct directions, i.e. up 
(+) or down (— ), at opposite boundaries, respectively. 
Hence, the PBC are released in that direction, while they 
are preserved in the remaining d — 1 directions. We 
calculated the GSs for the four possible combinations 

H — , — +, ++, . Three types of spin-regions can 

be distinguished. The first type can be flipped chang- 
ing a single boundary condition, i.e. from ++ — > + — . 
We call these spin region strong controllable. The sec- 
ond type is just called controllable, if it can be flipped 
by changing both boundary conditions. The third type 
forms fixed, stable islands, unaffected of any boundary- 
condition change. Examples for regions of such spins are 
shown in Fig. [T] 
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FIG. 1: Fixing the boundary spins (black bars on the left and 
right of every sketch), forces controllable spins to the direction 
of the boundary-sign. Areas with Si — +1 are encoded in gray, 
areas with spins Si = — 1 are white. The d — 1-dimensional 
surface between spins of opposite direction is called domain 
wall (black lines). In some regions, the random fields may 
freeze the spins into boundary-independent stable islands. Is- 
lands may include opposite-directed islands, etc. Hence, is- 
lands may also interfere with the domain wall. The dashed 
lines signal the border of islands within an area of equally- 
oriented spins, and, at the top and the bottom, the periodic 
boundary conditions in the remaining d — 1 directions. 



We also compare two kinds of droplet excitations of the 
GS. In both cases, the starting point is the GS calculated 
for full PBC of a realization. 

• The first type of excitation is obtained by fixing 
(L/3) d spins in the cneter opposite to the their 
ground-state orientation, respectively. This is in- 
spired by the approach of Ref. for a disordered 
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solid-on-solid model. This effect can be achieved 
conveniently by applying strong local fields in 
the desired direction for the fixed spins. Also, we 
fixed the spins on hyper-planes of the boundary 
in parallel to the GS orientation. 20 Note that the 
spins inside the bulk area create a contribution to 
the droplet energy via the local fields. To exclude 
this unwanted effect, we set the local fields e in the 
bulk region to zero for the intitial GS calculation. 
Via a recalculation of the GS of the modified sys- 
tem, this leads to a large excitation with respect to 
the original realization, not including the bound- 
ary, i.e., it is impossible that the fully system flips 
over. Below, we refer to these excitations as bulk- 
induced droplets, see Fig. [2j From the definition it 
is clear that these excitations involve 0(L 3 ) spins. 
Therefore, they come very close to the definition of 
droplets in the droplet theor y 20 ' 21 . To the knowl- 
edge of the authors, such droplets have not been 
studied for the RFIM so far. 

• The other type of excitation, called single-spin- 
induced droplets, consists in flipping only the 
very centered spin and freeze it anti-parallel to 
its ground-state orientation, again including fix- 
ing boundary spins in parallel to the GS. In the 
same way as for the bulk-induced droplets, this 
is achieved by applying strong local fields. The 
droplet created in such a way will include the cen- 
ter spin, but not the boundary. Theses singl-spin- 
induced droplets will be usually smaller than the 
bulk-induced droplets. Such excitations have been 
studied in d = 3 so far^°- 

Hence, for each realization {hei] of the disorder, 
we obtained seven different (ground-state) configura- 
tions for diffcrt types of boundary conditions/constraints 

(PBC, (++), ( ), (+-), (-+), bulk-induced droplets 

and single-spin-induced droplets). Technically, to obtain 
the droplets, we extracted the absolute differences be- 
tween the spin configuration of two or more GSs via linear 
combinations of the configurations. Spins with the same 
properties add to the same value and form connected 
clusters. These clusters were obtaineded using a breadth- 
fist-search algorithm. In this way, the calculation of ge- 
ometric properties of the clusters was very convenient. 
For the droplets, we measured the domain wall enclosing 
the droplet, i.e., the droplet surface. In the case of dif- 
ferent boundary conditions, we are interested in different 
definitions of the domain walls, i.e. separating uncon- 
trollable spins vs. controllable spins, strong-controllable 
spins vs. controllable spins and islands, allowing us to 
measure different related fractal exponents: 

The surfaces are usually not flat or smooth, instead a 
any type of disorder- averaged surface A exhibits asymp- 
totically a fractal scaling behavior of the form 




forced spins 



FIG. 2: Spin-wise difference of the pure GS and a droplet is 
shown. Unchanged spins are shown as white areas (pi — Sj = 
0), changed spins are displayed in gray (|sj — s'i\ = 2). 
The excitation is generated, freezing the spins on the bound- 
aries to their pure GS orientation, while the spins on a re- 
gion in the center are frozen inverted to their ground-state 
orientation (black dot for the single-spin-induced droplets or 
dashed square for the bulk-induced droplets). Under these 
frozen-spin constraints a new GS is calculated. This new con- 
figuration is an excitation with respect to the original GS. 
The resulting droplet may spread to an arbitrary shape with 
a fractal surface. It also may contain islands. 



d being the corresponding fractal dimension and c a 
constant. Dependent on the different boundary and 
droplet conditions several fractal exponents can be de- 
duced. Partially following the definitions of Ref. Q we 
measured 

d a : Surface exponent as dimensionality of the (hyper) 
surface of unchanged spins within two GSs of H — h 
and H — boundary conditions. Technically, we cal- 
culate Xi = \sf + — | = {0,2} and count the 
bonds between the 0-cluster (unchanged spins) and 
the 2-cluster (changed spins). 

di~. Incongruent interface exponent as dimensional- 
ity of the surface of unchanged spins within the 
three GS of ++, +— and boundary condi- 
tions. This domain-wall does not include parts 
of any stable islands. In detail, we calculate 
x, = s" + sf + - 2sr+ + 4 = {2, 4, 6}. It results 
to 2 if Si is strong-controllable. Spins with xi — 4 
belong to a stable island and Xi = 6, if Si is just 
controllable. For the incongruent boundary the 
bonds only between the 2-cluster and 6-cluster are 
counted. 

dj: Exchange stiffness exponent, defined as the singed 
sum of broken bonds, counted positive in the H — 
and — h configuration and negative for H — h and . 



A o =c L d % (3) 



dB- fractal exponent of the surface between all flipped 
and unflipped spins , i.e., number of bonds between 
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these, for bulk-induced droplets, including islands. 
d^: same as above, but excluding islands 

d\: fractal exponent of the surface between all flipped 
and unflipped spins, induced by a single-spin- 
undiced excitations, including islands. 

di : same as d±, but excluding islands. 

In the case of droplet excitations, we measured addi- 
tionally the disorder-averaged volume V, i.e., the average 
number of spins in the cluster of droplet spins 

Additionally to the geometric properties of domain 
walls and droplets we are interested in the stiffness ex- 
ponent 9, comparing three different types of excitations. 
The first approach is based on the symmetrized stiffness 
£ defined by Middleton and Fisher^ i.e., the disorder- 
averaged symmetrized sum of the boundary-condition de- 
pendent energies. 

£ = + E-+ - E ++ - E — ) /2 , (4) 

where E pq is the GS energy for boundary condition pq G 

{++, , H — , — h}, and (•) denotes the disorder average. 

A detailed picture of resulting configurations can be 
found in Fig. [TJ Close to criticality, the average stiffness 
can be assumed to scale as£ 

£(L)~L*. (5) 

9 denotes the stiffness exponent. According the droplet 
theory, the energy for other types of excitations of or- 
der of system sizes should scale with the same exponent. 
Hence, for the bulk induced droplets, one should be able 
to observe 

AE B (i)~i 9 , (6) 

where AEb is the average of the excitation energy of the 
droplet, i.e., the energy of the droplet configuration for 
the orginal value of the fields minus the GS energy for 
the same original values of the fields. 

Note that the single-spin-induced droplets tend to be 
small, i.e., they are not of order of system size. Hence, 
one cannot directly measure 9 from the scaling of the 
droplet energy. Instead, for the third approach, we fol- 
lowe the arguments of Ref. 3C|. Therein, it is shown 
that the distribution of single-spin-induccd droplet radii 
scales with p(R) ~ R . 

IV. RESULTS 

We performed exact ground-state calculations for 
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TABLE I: Fit parameters for the scaling of the different sur- 
face definitions, i.e. pure surface, £j and the incongruent 
parts of the domain walls. 



realizations of disorder, where the largest size exhibit 
the smallest number of realizations, respectively. Note 
that for d = 4, we studied L = 45 only for the bulk- 
induced droplets, while we did not included results for 
bulk-induced droplets in d > 5. In general, due to com- 
puter main memory restrictions we are limited to system 
sizes below 5 million spins. We start our analysis with 
the geometric properties of the domain walls. 

A. Domain walls 

The surfaces of the three different defined boundary- 
induced domain walls scale with plain and very clear 
power laws. Exemplary, the scaling of the simple domain- 
wall, i.e. (++)/(H — ) (yielding the fractal exponent d s ) 
is shown in Fig. [3] Error bars^i where obtained as stan- 
dard error bars from the empirical variance. The other 
plots look quite the same with the same precision. The 
scaling exponents d s , dj and di were found with high 
statistical accuracy. They are stated in Tab. HI Note that 
the upper limit for any fractal dimension is d, hence the 
result for d s at d = 7 is an artifact created by the small 
range of sizes which is accessible at this high dimension. 
Also, the obtained exponents depend on the fit range, 
reflecting possible systematic correction to the limiting 
scaling behavior of Eq. ©. To estimate such system- 
atic errors, we have performed fits for different ranges, 
leading to the final results which are displayed, c.f. Tab. 
IVII1 Given this accuracy, the results for dj may dj agree 
or differ, in particular in large dimensions, see discussion 
in Sec. El but d s differs, comparable to the previously 
obtained results^i 1 - for d = 3, 4. 



We continue calculating the stiffness exponents. First, 
the well established ansatz according to Eq. is used 
for d = 3,4, 5, 6, 7. Our data show well behaved power 
laws in each dimension, see Fig. [4] The fits of the data 
points follow Eq. ([5]), the parameters are stated in Tab. 
HI! in particular the resulting estimate for the stiffness 
exponentn 9 is denoted as 9& w . Varying the ranges of 
fitted sizes leads to the final estimates, again stated in 
Tab.rVHl 
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FIG. 3: (color online) Simple surface of domain walls. Most 
error-bars are smaller than the symbol size. The lines show 
power-law fits according to Eq. ((3J- 
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TABLE II: Fit parameters for 6 according to Eq. ((5j 
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FIG. 4: (color online) Stiffness, as defined in Eq. @. Most 
error bars are smaller than the symbol size. Lines display 
results of fits to power-laws according to Eq. l[5|. 



B. Droplets 

Concerning the droplet behavior, we start with the 
fractal surface properties of the droplets. 

We start by using a scatter plot to obtain the rela- 



tion between surface and volume of single-spin-induced 
droplets. Theoretically, it should follow a power law of 
the form 



A ~ V dl/d . 



(7) 



The fractal surface exponent with and without stable is- 
lands can be obtained this way, depending on whether 
the islands are included in the calculation of A. And 
indeed, the results exhibit clear power laws, see Fig. [5]). 
Additionally, for each dimension the data points for other 
system sizes L scatter around the same line, respectively. 
The fit parameters are listed in Tab. IIII1 for fixed sys- 
tem sizes. The fractal exponents including or ignoring 
enclosed islands do not differ significantly. We have ob- 
served slight changes of the value of d-y when varying the 
system size L (unless using very small system sizes where 
the fitted value of d\ is much smaller). Hence, the final 
values we quote, c.f. Tab. I VIII are slightly different and 
involve larger error bars than the pure statistical error 
bars. When approaching large dimensions d s gets close 
to the dimension of the system, in particular right at the 
upper critical dimension. The result for the fractal di- 
mensions is, given the unknown correlations to scaling, 
in fair agreement with the results for d s as shown in Tab. 

in 
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FIG. 5: (color online) Scatter plot of the enclosing surface, 
i.e., not counting islands inside the droplet, of single-spin- 
induced droplets as function of their volume. The surface is 
scaled up with a factor of k — 1, 2, 4, 8, 16 for d — 3, 4, 5, 6, 7 
to separate the data visually. 



The bulk-induced droplets have the disadvantage that 
a smaller effective range of sizes if accessible: The min- 
imal system size to generate such a droplet excitation 
must be L — 6. In linear direction, one spin is used to fix 
the boundary and at least another one is needed at the 
very center to form the "large core". Though, only two 
spin in each direction are left to form the droplet. In the 
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TABLE III: Fit parameters for the scaling of the enclosing 
surface and ragged surface for single-spin-induced droplets, 
calculated in a large system size with high sample number, 
system sizes as in Fig. [5] 
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TABLE IV: Fit parameters for the scaling of the enclosing 
surface and ragged surface for bulk-induced droplets. For 
d > 4, data is not available in high quality for large-enough 
system sizes. 



same way, also for larger sizes, the volume accessible for 
the droplets to form is smaller, leading to stronger finite- 
size corrections, see below. This means, for d > 5 the 
range of system sizes (L > 6) is too small to observe the 
leading scaling behavior. Hence, we restrict our analysis 
for the bulk-induced droplets to d = 3, 4 and d — 5. 

In Fig.|6]a scatter-plot of the enclosing surface as func- 
tion of the volume can be seen. For the fits the ansatz 
according Eq. ([7]) was made, but using dB instead of d±. 
From the large number of data points the fit parameters 
can be obtained, by using the data points for all system 
sizs in one fit, with high statistical accuracy. The fit pa- 
rameters can be found in Tab. II VI Given the small ranges 
of system sizes here and unknown corrections to scaling, 
the agreement with the results for the single-spin-induced 
droplets is fair. 



Systematic finite-size corrections are negelected by this 
single fit of the scatter-plot data. Hence, to estimate 
these corrections, leading to the results stated in Tab. 
IVII1 we proceeded as follows: Instead of observing the 
scatter-plot, the mean of the enclosing surface can be 
analyzed as function of the system size. A power law 
as A B (L) = c B {L - L ) dB leads to the best fits. The 
length-scale correction Lq covers finite-size effects. Of 
course, an ansatz using a correction term CBL ds (l + 
C2U 1 ) can be used too, but the correction term involves 
one more parameter. Anyway, fitting Ab{L), leads to 
similar exponents compared to those as obtained from 
the scatter plot, see Tab. |V] Again, the results for the 
different approaches for surface measurement do not lead 
to significantly different results. 

Next, we discuss energetic properties of droplets, aim- 
ing at estimating the stiffness exponent 8. 




FIG. 6: (color online) Scatter plot of the enclosing surface 
of bulk-induced droplets as function of their volume for d — 
3, 4, 5. The lower right labels show system sizes of the 3d 
data, the upper left show those for 5d. 4d is unlabeled. 
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TABLE V: Fit parameters for the scaling of the enclosing 
surface and ragged surface for bulk-induced droplets when 
scaling the surface as a function of system size. For d = 6,7, 
a large-enough range of system sizes is not available. 



The singe-spin droplets are known to tend to be very 
small. 30 This property prohibits the direct measurement 
of the stiffness exponents from the relation of droplet en- 
ergy to linear droplet size. Nevertheless, the distribution 
of droplet radii R follows 30 approximately a power law 
P(R) ~ R~ e , when the data is binned logarithmically. 
Hence, we obtained the distribution of droplet radii. For 
d = 3, 4, 5 a logarithmic binning of the data points is 
possible. For larger dimensions, the achievable systems 
sizes seem to be too small to get a histogram of sufficient 
statistical quality. The histograms can be seen in Fig. 
[7] There exist regions in which the radii distribution fol- 
lows R~ e , when using the values #dw obtained from the 
domain-wall measurements. This confirms somehow the 
results of Ref<^ for higher dimensions d = 4,5. Neverthe- 
less, the distribution itself is not sufficient to determine 
the stiffness in a meaningful way if the value of 8 was not 
known from other sources. 



For the bulk-induced droplets, their length scale is 
fixed to be O(L), as assumed in the droplet theor y 20 ' 21 . 
Hence, one could hope that the scaling of the droplet en- 
ergy is governed by the (stiffness) exponent 8. We looked 
at the scaling of the difference between the GS energies 
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of the original system and the GS for the bulk-induced 
droplets (calculated with the original set of random fields, 
respectively). The resulting average bulk droplet energies 
AEb are shown in Fig. [5] for d = 3, 4. A clear curvature 
is visible, hence strong finite-size corrections to Eq. ([6]). 
This can be explained by the fact that, due to the exten- 
sive size of the bulk area, the effective volume which is 
accessible for the droplet to arrange is much smaller. In 
particular this means, that for higher dimension d > 4 
the linear sizes L which are accessible are to small to 
come even just near to the final scaling behavior. Hence, 
we have restricted ourself for the energetic properties of 
the bulk-induced droplets to dimensions d — 3 and 4. To 
include corrections to scaling, we fitted the data can by 
using 

AE B = a(L - Lo) 9b , (8) 

see Tab. IVII The minimum system size L m - ln included 
in the fit was chosen such that the fitting quality (mea- 
sured by the weighted sum of square residuals per degree 
of freedeom, denoted as WSSER/ndf in gunplot) was ac- 
ceptable, i.e., around 1. To furthermore estimate system- 
atic corrections, we also performed fits for larger values 
of ^min, respectively, leading to the final results for 9b as 
stated in Tab. EED Within error bars, the values for the 
droplet stiffness are compatible with the results obtained 
for the domain walls. Hence, it appears that indeed the 
basic assumptions of the droplet theory is true that dif- 
ferent types of excitations are universally described by 
the same exponents. 




R 



FIG. 7: (color online) Number of single-spin-induced droplets 
with a radius in a specific interval for d — 3, L = 96, d = 
4, L — 24 and d = 4, L — 13. The latter down scaled by a 
factor 10 and 20, respectively. The lines are power laws using 
the stiffness exponents. 



d Lmj n a Lq 8b 

3 15 5.34(15) 2.42(10) 1.501 (7) 

4 12 7.41(18) 2.75 (5) 1.91 (1) 

TABLE VI: Fit parameters for the scaling of the droplet en- 
ergy according to Eq. ([B} • The fits were performed for a range 
of system sizes rangeing from L m i n to the largest system sizes 
considered here, respectively. 




L 



FIG. 8: Scaling of the excitation energy AEb of bulk-induced 
droplets for d — 3 and d = 4. Lines show fits to power-laws 
AE B ~ (L - L ) s , with 9 = 1.50 (3d) and 6 = 1.91 (4d). 

C. Scaling of ground-state energy 

There seems to exist a third alternative way of de- 
termining the stiffness exponent, not related to exter- 
nally induced domain walls and droplets. For spin 
glasses in finite and low dimensions with PBC, it has 
been conjectured^ - and numerically observed for two 
dimensions^ that the finite-size behavior of the (total) 
GS energy El of the unperturbed system is given by 

E L = E oa + aL e . (9) 



The explanation is that the PBC induce "hidden" domain 
walls (with respect to free boundary conditions), which 
dominate the finite-size corrections above higher order 
contributions. Hence, we have reanalyzed the data of 
previous work ^ 10 ' 13 The results of 6 from fitting ([9]) to 
the GS energies can be found in Fig. |H1 its exponents are 
listed as $e in Tab. I VIII Comparing with the results for 
the stiffness exponent #dw, the assumptions seems to be 
indeed true for d < d u — 6, maybe also at d = d u (given 
that we state here only statistical error bars, hence true 
error bars are likely larger). The assumption of Eq. © 
is certainly not true above the upper critical dimension. 
It seems rather that the exponent for GS energy finite- 
size correction attains some mean-field value for d > d u . 
Interestingly, #e agrees within error bars with "f /v in all 
dimensions d = 3, . . . , 7. 
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FIG. 9: (color online) Scaling of the ground-state energy at 
the critical point. The lines show fits to Eq. ([9|. 



V. CONCLUSION AND DISCUSSION 

To calculate ground states of the RFIM numerically, 
we applied a well-known mapping to the maximum- 
flow problem. Using efficient polynomial-time-running 
maximum- now/minimum-cut algorithms, we were able 
to study large systems sizes up to 5 x 10 6 spins in exact 
equilibrium at T = 0. 

Comparing the GS configurations obtained from dif- 
ferent constraints leads to suitably defined domain-wall 
and droplet excitations. We obtained fractal surface ex- 
ponents and the stiffness exponent in d = 3,4,5,6,7 
for these excitations, see final results in Tab. IVIII The 
values for domain walls stated in the literature 7 -^ for 
d = 3,4 could be recovered. For all dimensions, differ- 
ent types of excitations are described, within error bars, 
by the same value of the stiffness exponents, as assumed 
by the droplet theor y 20 ! 21 In particular, (bulk) droplet 
excitations have been calculated for the RFIM for the 
first time and found to be compatible with domain walls. 
Also, the result for 6*d w is according a scaling relation 7 
equal to dj — which is fulfilled within error bars for 
all dimensions. Furthermore, the scaling inequalit y 44 i 45 
d/2- fi/v <6 < 6/2 is fullfilled in all dimensions. 

On the other hand, the hyperscaling relation ([1]) is ful- 
filled within error bars only for d = 3,4,5, but not for 
d = 6, 7. Usually, hyper-scaling relations, involving the 
dimension d, are not valid im the mean- field regime, i.e., 
above the upper critical dimension. Nevertheless, moti- 
vated by our results for d = 7, it is compatible with the 
data if we change the hyperscaling relation in the fol- 
lowing way: We know from earlier work in high dimen- 
sions that above the upper critical dimension d u when 
performing a data collapse, a rescaling of the linear sys- 
tem size according to L — > L d / du has to be performed. 46 
Since the length rescaling enters into scaling relations 
via the correlation lengtbJ 7 -^ £l/L ~ l(L x / u (h - h c )), 
the correlation length exponent v needs to be replaced 



h c 

P 
7 

Q 
V 



dj-l/v 



d s 
di 
d B 
dj 
di 



d- 



d- 



d = 5 d — 6 d 



2.28(1) 
0.017(5) 
1.98(7) 


1.37(9) 
1-4(2) 



4.18(1) 
0.13(5) 
1.57(10) 


0.78(10) 
2.0(4) 



6.02(2) 7.78(1) 
0.25(1) 0.50(5) 
1.3(1) 1.07(5) 

0.60(3) 0.50(5) 
2.3(3) 2.1(5) 



9.48(3) 
0.50(5) 
1.0(2) 


0.47(5) 
2.2(5) 



1.49(3) 
1.44(2) 
1.4(1) 
1.51(2) 



1.81(6) 
1.75(2) 
1.65(15) 
1.8(1) 



2.03(2) 2.42(2) 
2.15(1) 2.60(2) 
2.1(1) 2.5(2) 



2.42(2) 

3.5(1) 

3.5(1) 



2.37(2) 3.92(6) 4.9(2) 5.8(8) 6.9(2) 

2.32(1) 3.73(10) 4.9(1) 5.9(1) 6.9(1) 

2.35(2) 3.79(4) 4.6(1) 

2.14(3) 2.93(2) 3.79(2) 4.54(2) 5.5(1) 

2.25(2) 3.23(1) 3.92(1) 4.5(2) 5.4(2) 



TABLE VII: Previous results 7 - ' 11 ' 13 above the big line and 
new, final results below. The upper half of the final results 
contains different estimates for the stiffness exponent, while 
the bottom part contains measured for the fractal proper- 
ties of domain walls. Note that the values stated here differ 
usually from the values given in the lists of fit parameters 
since they the former include systematical errors, which were 
estimated by varying the range of fitted system sizes and ob- 
serving the change of the resulting parameters. 



with v* = vyi^du/d when performing a data collaps. If 
we transfer this change to the hyper-scaling relation, it 
reads 9 = d — (2 — a)/v* . Inserting in d — 7, d u = 6, 
a = and i^mf = 0.5, this leads to 6^ = 3.5 with- 
out any further assumptions. This matches very well to 
#( 7 ) = 3.5(1) as found in this study. On the other hand, a 
simple rescaling in d — 6 appears not possible to us, since 
a replacement L — V L ln^L was found previously^ so 
rescaled exponents can not be obtained in a simple way, 
hence a different kind of hyperscaling seems to be neces- 
sary here. 

Also the fractal properties of the droplets and domain 
walls found for lower dimensions extent to larger dimen- 
sions. In particular, the "geometrically" measured val- 
ues d s , d\ and g?b agree for all dimensions within error 
bars. Finally, the values for dj and dj, measuring differ- 
ent fractal properties of the domain walls, may or may 
not be equivalent, as already discussed in previouslyii 1 - 
for the cases d = 3 and d = 4, were also a possible scal- 
ing relation dj = dj + /3/u was presented. For d > 6 this 
relation is clearly not fullfiled here, which may be also 
due to the failing of hyperscaling at and above the upper 
critical dimension. 
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